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Main  Scientific  Objective 


•Inversion  of  the  bottom  parameters  through  reflection 
phase  shift  of  sea  bottom 

Step  1 :  Extract  modal  depth  function  from  acoustic  data 

Step2:  Estimate  km  from  extracted  modal  depth  function 

Step3:Calculate  reflection  phase  shift  of  sea  bottom 

Step4:  Estimate  bottom  parameters  by  use  of  reflection 
phase  shift  of  sea  bottom 

•Study  the  possibility  of  extracting  modal  depth  function 
from  reverberation  data  collected  by  a  VLA 


Relation  between  the  phase  shift  of 
sea  bottom  and  modal  wavenumber 


The  bottom  reflection  phase  shift  can  be  retrieved  if 
the  real  part  of  the  modal  wavenumber  and  the  SSP 
in  water  column  are  known. 
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The  relation  used  for  extraction  of  modal  depth 
function 


•The  reverberation  field  received  by  a  hydrophone 


•The  expression  of  the  scattering  field  received  by  a  VLA 


where 


The  relation  used  for  extraction  of  modal  depth 
function 


CSDM  of  reverberation 
VLA  data 


SVD  of  the  CSDM  of 
reverberation  VLA  data 


*  These  relations  are  going  to  be  verified  in  detail  later 
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Estimation  from  extracted  modal  depth  function 
(shooting) 


Assume  a  k 


Measured  SSP 


Finite 
difference 
equation 
over  depth 


odal  depth  function 


Extracted 
modal  depth 
function 


F(km)  =  £[$,  .(z,)-^(z,)]: 
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De-noised  modal  depth  function 


Experiment  setup 
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lKg  explosive 
source 


•  Not  used  in  this  paper 
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Extracting  modal  depth  function  from  reverberation 
data  method  I:  Different  Time-segment  Method 
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Shooting 


De-noised  <I> 


Extracting  modal  depth  function  from  reverberation 
data  method  II:  Different  Frequency  Method 


fj=fO+[-NAf,(-N+l)  Af, . . .,N  Af| 
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Shootin 


De-noised  <3>m 
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^  |  Example  of  extracted  mode  depth  function 
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*  ^Different  time- segment  method 


•f0=100Hz 


Relative  singular  value  Extracted  modal  depth  function 
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*  extracted 
—  shooting 


Example  of  extracted  mode  depth  function 

•Different  Frequency  method 
•f0=200Hz,  Af=lHz,  N=20 


Relative  singular  value  Extracted  modal  depth  function 
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^  |  Retrieved  modal  phase  and  verify 


Example 

m 

k 

m 

cf> 

Paw  i 

m 

P 

100Hz 

Multi-r 

C 

4 

0.3980 

-0.072 

3.0698 

0.2754 

11.016 

5 

0.3896 

0.632 

3.774 

0.3436 

10.984 

200Hz 

Multi-f 

6 

0.8066 

-0.689 

2.453 

0.2238 

10.961 

9 

0.7833 

0.581 

3.723 

0.3276 

11.36 

10 

0.7734 

1.11 

4.2516 

0.3630 

11.71 

From  WKB  approximation,  the  reflection  phase  shift  is  linear  with 
grazing  angle  when  the  grazing  angle  is  far  smaller  than  1. 


arg[fl(0 )]  = 


where 


Inversed  geoacoutic  parameters  of  sea  bottom 


What  does  P  value  means? 


Table  2  Hamilton’s  typical  sediment  parameters  and  P[9] 


Type 

Sediment  parameters 

P 

Q 

Critical 

Density(g 

/cm3) 

Sound 

speed(m/s) 

K 

angle(°) 

1 

2.03 

1836 

0.479 

7.3 

0.22 

34 

2 

1.97 

1753 

0.510 

7.9 

0.34 

29 

3 

1.87 

1697 

0.670 

8.5 

0.63 

26 
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0.692 
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0.80 

24 
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0.90 
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1623 

0.673 
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20 
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1.58 
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13 

0.52 
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1.47 
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0.090 

20 

2.27 
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1520 
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A  specific  P  value  corresponds  a  specific  type  of  sediment! 
The  inverted  P  value  is  about  1 1 ,  it  is  very  close  to  the  P  of 
type  6(sand  silt).  So  we  know  the  density  and  sound  speed 
are  close  to  1.7  and  1623,  respectively. 


Error  Analysis 


where 


A (m)  ~Sm-Ak 
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Summary 


The  modal  depth  function  can  be  extracted  by  use  of  reverberation 
VLA  data. 

The  modal  wavenumber  can  be  estimated  from  extracted  modal  depth 
function  by  shooting  method.  The  error  of  estimated  modal  wavenumber 
is  about  10'4~10-3. 

A  consistent  P  value  can  be  obtained  from  the  retrieved  modal  reflection 
phase  shift. 

From  the  inverted  P  value,  the  sediment  type  of  the  sea  bottom  where  the 
reverberation  experiment  was  conducted  in  ASIAEX  is  corresponded  to  the 
‘type  6’  of  the  Hamilton’s  classification  (sand  silt:density=1.7g/cm3,  sound 
speed=1623m/s). 

In  this  paper,  we  only  extracted  higher  order  modes,  we  think  it  is 
because  the  VLA  is  short  and  span  only  the  upper  half  of  the  water.  of  0c%r 
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Ambiguity  of  inversion  by  using  matched  field 
process 
Example: 


VLA  :  l~103m  with  spacing  lm 
Sea  Depth:  105m 
Sea  bottom:  Half-infinite 
Parameters  for  Measured  field: 
density = 1 .7  g/  cm3 
Sound  speed=1623m/s 
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Ambiguity  of  inversion  by  using  matched  field 
process 


Ambiguity  surface  for  the 
correlation  coefficient 


Comparing  with  our  result,  MFP  inversion  has  a  large 
ambiguity.  In  our  method,  the  ambiguity  is 
significantly  reduced  by  taking  Hamilton’ s  emperical  j?- of  °c<>\ 
data  as  a  “constraint”. 
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